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Abstract. We consider the 3D quantum many-body dynamics describing a dilute bose 
gas with strong confining in one direction. We study the corresponding BBGKY hierarchy 
which contains a diverging coefficient as the strength of the confining potential tends to 
oo. We find that this diverging coefficient is counterbalanced by the limiting structure of 
the density matrices and establish the convergence of the BBGKY hierarchy. Moreover, we 
prove that the limit is fully described by a 2D cubic NLS and obtain the exact 3D to 2D 
coupling constant. 
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1. Introduction 
It is widely believed that the cubic nonlinear Schrodinger equation (NLS) 

idtcj) = L0+ in 

where L is the Laplacian —A or the Hermite operator —A + u"^ \xf , describes the physical 
phenomenon of Bose-Einstein condensation (BEC). This belief is one of the main motivations 
for studying the cubic NLS. BEC is the phenomenon that particles of integer spin (bosons) 
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occupy a macroscopic quantum state. This unusual state of matter was first predicted 
theoretically by Einstein for non-interacting particles. The first experimental observation of 
BEC in an interacting atomic gas did not occur until 1995 using laser cooling techniques 
[HEQ]. E. a. Cornell, W. Ketterle, and C. E. Wieman were awarded the 2001 Nobel Prize in 
physics for observing BEC. Many similar successful experiments [HI |36l [52] were performed 
later. 

Let t G M be the time variable and = (ri,r2, ■■■,r]\f) G M."'^ be the position vector of N 
particles in M". Then BEC naively means that the A^-body wave function ip^it, r^) satisfies 

N 

up to a phase factor solely depending on t, for some one particle state ip. In other words, 
every particle is in the same quantum state. Equivalently, there is the Penrose-Onsager 
formulation ^46j of BEC: if we define 7^'* to be the A;-particle marginal densities associated 
with by 

(1.1) 75J'^(t, Tk, r'fc) = j ij^{t, Ffc, rjv_fc)V^7v(^, r'^, rN-k)drN-k, r^, 4 G M"'' 

then, equivalently, BEC means 

k 

(1-2) jP{t,r,-r',)^l[^{t,rjMt,r'^). 

i=i 

Gross [331 131] and Pitaevskii [17] proposed that the many-body effect should be model by 
a strong on-site interaction and hence the one-particle state if should be modeled by the a 
cubic NLS. In a series of works [HUlllEIlISSllMlESlEniEZllIIlIIZ] , it has been proven 



rigorously that, under suitable assumptions on the interaction potential, relation (1.2) holds 
in 3D and the one-particle state cp satisfies the 3D cubic NLS. 

It is then natural to believe that the 2D cubic NLS describes the 2D BEC as well. However, 
there is no BEC in 2D unless the temperature is absolute zero (see p. 69 of [13] and the 
references within). In other words, 2D BEC is physically impossible due to the third law of 
thermodynamics. In a physically realistic setting, 2D NLS can only arise from a 3D BEC 
with strong confining in one direction (which we take to be the z-direction). Such an effective 
3D to 2D phenomenon has been experimentally observed [2HI |53l [191 ISHl [18]. (See [6] for 
a review.) It is then natural to consider the derivation of the 2D NLS from a 3D A^-body 
quantum dynamic. Combining [H [21 |T7] suggests a route of getting the 2D NLS from 3D. 
First, a special case of Theorem 2 in [T7] establishes the 3D cubic NLS 

(1.3) idt'p = -A^<p+ {-dj + uj^z'^)<p+\<pf<p, {x,z)eR^^^ 

from the 3D A^-body quantum dynamic as a A^ — ?■ 00 limit. Then the result in [H [2] shows 



that the 2D cubic NLS arises from equation (1.3) as a a; —;■ 00 limit. This path corresponds to 



the iterated limit (lim^_5.oo liniAr_>.oo) of the A^-body dynamic, thus the 2D cubic NLS coming 
from such a path approximates the 3D A^-body dynamic when u is large and A^ is infinity. 
In experiments, it is fully possible to have A^ and co comparable to each other. In fact, A^ 
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is about 10^ and u is about 10^ in [28], |53l [351 IIH]- In this paper, we derive rigorously the 
2D cubic NLS as the double limit (limAr^^_>oo) of a 3D quantum iV-body dynamic directly, 
without passing through any 3D cubic NLS. It is elementary mathematical analysis that 
lim^_>oo liniAT-s-oo and \mi]s[^u)^oo are topologically different and one does not imply each other. 
Let us adopt the notation 

and investigate the procedure of laboratory experiments of BEC according to [28], 153] [191 ES] 

[ig. 

Step A. Confine a large number of bosons inside a trap with strong confining in the z- 
direction. Cool it down so that the many-body system reaches its ground state. It is 
expected that this ground state is a BEC state / factorized state. To formulate the problem 
mathematically, we use the quadratic potential |-| to represent the trap and 

to represent the interaction potential. We use the quadratic potential to represent the 
trap because this simplified yet reasonably general model is expected to capture the salient 
features of the actual trap: on the one hand the quadratic potential varies slowly, on the 
other hand it tends to cxo as |x| — )■ oo. In the physics literature, Lieb, Seiringer and Yngvason 
remarked in [H] that the confining potential is typically ~ \xf in the available experiments. 
The review [6] on [28] [53] [191 [SH] [18] also mentioned that the trap is harmonic. We use Va (r) 
to represent the interaction potential to match the Gross-Pitaevskii description [33l [311 HZ] 
that the many-body effect should be modeled by an on-site self interaction because Va is 
an approximation of the identity as a — )■ 0. This step then corresponds to the following 
mathematical problem: 

Problem 1. Show that, for large N and large u ^ uoq, the ground state of the N-body 
Hamiltonian 

j=l l^i<j^N ^ ^ 

is a factorized state under proper assumptions on a and V. 

Step B. Switch the trap in order to enable measurement or direct observation. It is assumed 
that such a shift of the confining potential is instant and does not destroy the BEC obtained 
from Step A. To be more precise about the word "switch": in [THl [IB], the trap in the 
x-spatial directions are tuned very loose to generate a 2D Bose gas. For mathematical 
convenience, we can assume Uq becomes 0. The system is then time dependent. Therefore, 
the factorized structure obtained in Step A must be preserved in time for the observation of 
BEC. Mathematically, this step stands for the following problem. 
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Problem 2. Take the BEC state obtained in Step A as initial datum, show that, for large 
N and u, the solution to the N—body Schrddinger equation 

(1.5) 0.^„,„ = f (-iA,, + ^.|)v.„,„+ E J^v('-l-p)i, 

j=l ^ ^ l^i<j^N ^ ' 

is a BEC state / factorized state under the same assumptions of the interaction potential V 
in Problem{l\ 

We first remark that neither of the problems listed above admits a factorized state solution. 
It is also unrealistic to solve the equations in Problems [T] and [2] for large N. Moreover, both 
problems are linear so that it is not clear how the 2D cubic NLS arises from either problem. 
Therefore, in order to justify the statement that the 2D cubic NLS depicts the 3D to 2D 
BEC, we have to show mathematically that, in an appropriate sense, for some 3D one particle 
state ip fully described by the 2D cubic NLS 

k 

7!J,l(^,rfc;rfc) ~ JJ rj)^(t, r^-) as A^,u; ^ oo 
i=i 

(k) 

where "Jnuj are the fc-marginal densities associated with '?/'7v,a;- 

For Problem [l] (Step A), a satisfying answer has been found by Schnee and Yngvason. 
Let scat(W^) denote the 3D scattering length of the potential W. By pU Lemma A.l], for 
< /3 < 1 and a ^ 1, we have 



f ^ \a V ifO</3<l 

^ ^ [ascat(l^) if /3 = 1 

Consider (p^^ ^g, the minimizer to the 2D NLS energy functional 

(1.6) E^o,Ng= f {\V<j){x)\^ + ujl\xf\<j){x)\^ + 47rNg\<f){x)\^)dx 



subject to the constraint ||0||^2(]j2) = 1. The existence of this nonlinear ground state stems 
from the presence of the confining potential otherwise the nonlinear term is defocusing 

(as it is called in the NLS literature). 

Given parameters LOo,u,N,a, Schnee- Yngvason [19] define g = g{uJo,uj, N,a) and p = 
p{ijjo,u,N,a) by the two simultaneous equations (see (1.15) and (1.18) in [i9] ) 

-1 

-p = N 



def 
9 = 



They argue that this definition for g makes the 2D NLS Hamiltonian (1.6) relevant to the 
analysis of the limiting behavior of the ground state of (1.4) describing a dilute interacting 
Bose gas in a 3D trap that is strongly confining in the z-direction. (See also |54j for the case 
with rotation) 

The Gross-Pitaevskii limit means Ng ~ 1. We have liberty to fix the value of uq by 
scaling, so we take cjq = 1. Then the minimizer 4>u)o,Ng is fixed and hence p ~ iV. 
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In this paper, we consider Problem [2] (Step B) and offer a rigorous derivation of the 2D 
cubic NLS from the 3D quantum many-body dynamic. For the scahng of the interaction 
potential, we consider the case (called Region I in [19]) in which the term {^/ua)~^ dominates 
in the definition of g. Then 

1 ~ Ng ~ Na\/uj a ~ — 

This then implies that 

> log — ~ log 



so that our assumption that the term (a/cJo)"^ dominates in the definition of g is self- 
consistent. 

We will take for mathematical convenience a = {N^/uJ)~^ for Problem [2] (Step B). The 
Hamiltonian (1.4) then becomes 

N 

(1.7) ^^. = E(-^o+-%') + ]^ E {NV^rv{{NV^) 

j=l ^ l^i<j^N 

2/9 , 1 , T . , 1 1.1 1 , , . n , ■ n no . 2 



Let h{^z) = TT ^ so that h is the normalized ground state eigenfunction of —d1 + z 



1 



i.e. it solves (—1 — 9^ + z^^h = 0. Then the normalized ground state eigenfunction hi^{z) 
of — 9^ + u'^z'^ is given by huj{z) = uj^/^h{ijj^/'^z), i.e. it solves (— — dl + u'^z'^)h^ = 0. In 
particular, hi = h. 

We consider initial data that is asymptotically (as — 00, u — )■ 00) factorized in the 
z-direction and in the ground state in the 2;-direction; in particular we could take 

N 

^nA^^^n) = YlM^j)hMj) > ll0olU2(R2) = 1. 

Let 

(1-8) ^^,Jt,-) = e^*''"-^^,jO,-) 



denote the evolution of this initial data according to the Hamiltonian (1.7). We prove that 
in a certain sense, as — )■ 00, w — )■ 00, 

N 

(1.9) i>N,ujit^ ^n) ~ Yl ^^^^ ^i)^'^(^i) 

where 0(t) solves a 2D cubic NLS with initial data (poix). To make this statement more 
precise, we introduce the rescaled solution 

and the rescaled Hamiltonian 

^ 1 
(1-11) ^7V,. = E(-^-'.+^(-^i+^|)) + ]^ E ^NAr^-r,) 

j=l i<i<j<N 
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/3 {N^) 



and hence when ipiy^^{t) is given by (1.8) and ipj^^^ is defined by (1.10), we have 



The informal statement of convergence given by (1.9) becomes the informal statement 

AT 

^(^.rjv) ~ Jj0(t,a;j)/i(zj) 



(1.13) 



where (j){t) solves 2D NLS with initial data (p^ix). In fact, the convergence we prove is stated 
in terms of the associated density operators with kernels 



(1.14) 



The version of (1.13) that we prove is the convergence 

k 

7!J,l(t,r.,4) ^ l[<p{x,)h{z,)<P{x'^)h{z'^) 

in trace class, for each k > 0. 
We define 



:i.i5) 



v{(3) = max 



1-/3 iP-h ll^+l 1+i 
2/3 ' 1 - 1/3 ' 1-/3 ' 1 



2/3 



(see Fig. Q 

Our main theorem is the following: 



Theorem 1.1 (main theorem). Assume the pair interaction V is a nonnegative Schwartz 
class function. Let {7^'*^(t, r^; r'^,) } be the family of marginal densities associated with the 3D 
rescaled Hamiltonian evolution ip^^^^it) = c^^^'^''^ip^^^{fd) for some (3 G (0,2/5), (see (1.1), 
(1.11), (1.14)j. Suppose the initial datum ip ^{fS) satisfies the following: 

(a) i^N^^iGi) is normalized, that is, ||'0iv,a;(O)||L2 = 1, 

(b) V'Araj(O) asymptotically factorized in the sense that 



lim Tr 



7!i,l(0, xi, zi] x\,z[) - (f)Q{xi)(f)Q{x[)h{zi)h{z[] 



0, 



for some one particle state 0q G (M^) , 

(c) Away from the z-directional ground state energy, i'N^^iO) has finite energy per particle: 

sup JO), {Hn,. - Nu)i,^jO)) ^ C, 

uj,N 
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Then Vfc ^ l,t ^ 0, and e > 0, we have the convergence in trace norm (propagation of 
chaos) that 



lim Tr 

N,uj—^c<) 



iN^Jt^ Xfc> Zfc; x'^, Zfc) - Y[ (pit, Xj)(j){t, Xj)hi{zj)hi{zj] 



0, 



where v{j3) is given by (1.15) and (f){t,x) solves the 2D cubic NLS with coupling constant 
bo (^f \hi(z)\ dz) that is 



idt(t)=-l\A + h[ \h^{z)\Uz]\(t)\^(t) 



(1.16) 



with initial condition (0, x) = (f)o{x) and bo = j V (r) dr. 



in 



b2+1 



Theorem 1.1 is equivalent to the following theorem. 

Theorem 1.2 (main theorem). Assume the pair interaction V is a nonnegative Schwartz 
class function. Let {7^'|^(t, r^; r'^,) } be the family of marginal densities associated with the 3D 
rescaled Hamiltonian evolution i^j^^it) = e^^^^-'^ip ^ ^{Q) for some (3 G (0,2/5), (see (1.1), 
(1.11), ( 1.14[ ) ). Suppose the initial datum V'iv,a;(0) ^■^ normalized, asymptotically factorized 
and satisfies the energy condition that 
ic') there is a C > such that 



(1.17) 



Then ^ l,t ^ 0, and e > we have the convergence in trace norm (propagation of 
chaos) that 



lim Tr 



7!J,l(^, Xfc, Zfc; x'^, z'^) - JJ 0(t, Xj)(p{t, x'j)hi{zj)hi{zj] 



0, 



where v{/3) is given by (1.15) and (f){t,x) solves the 2D cubic NLS (1.16). 



We remark that assumptions (a), (b), and (c) in Theorem |1.1| are reasonable assumptions 
on the initial datum coming from Step A. In fact, if we assume further that 0o 



mmimizes 



the 2D Gross-Pitaevskii functional (1.6), then (a), (b) and (c) are the conclusion of 
Theorem 1.1, 1.3]. The limit in Theorem 11. ll which is taken as A^, a; — )■ oo within the 
subregion N ^ uj^(^)+'^ is optimal in the sense that if iV ^ a;2?~2^ then the limit of Vn,uj 
defined by (1.12) is not a delta function. 



The equivalence of Theorems |1.1| and |1.2| for asymptotically factorized initial data is well- 
known. In the main part of this paper, we prove Theorem|1.2|in full detail. For completeness. 



we discuss briefly how to deduce Theorem lA from Theorem L2 in Appendix [Bj 

The main tool used to prove Theorem L2 is the analysis of the BBGKY hierarchy of 

a.s N,u — 7- oo. With our definition, the sequence of the marginal densities 
associated with ipj^^^ satisfies the BBGKY hierarchy 



~{k) 



N 



k=l 

N 



k=l 
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Figure 1. A graph of the various rational functions of /3 appearing in (1.15). 
In Theorems 1.1, 1.2, the hmit {N,u) — )■ oo is taken with N > As 
shown here, there are values of (3 for which v{(3) ~ 1, which allows iV ~ w, 
as in the experimental paper [2S1 ESI IB]- We conjecture that Theorems 
hold with (1.15) replaced by the weaker constraint v{/3) = for all 



1.1 



1.2 

^ 1. 



(1.18) 



N -k 



+ E 



rk+l 



-dl + z 



1 \Vn,u^ {Tj 



2 

'j ' 'N,u. 



+ 



N 



E 

Kj 



rk+l] 



Vr 



~{k) 

> 7iV,a; 



In the classical setting, deriving mean-field type equations by studying the limit of the 
BBGKY hierarchy was proposed by Kac and demonstrated by Landford's work [H] on 
the Boltzmann equation. In the quantum setting, the usage of the BBGKY hierarchy was 
suggested by Spohn [SI] and has been proven to be successful by Elgart, Erdos, Schlein, and 
Yau in their fundamental papers [231 [211 12S1 [221 12Z] which rigorously derives the 3D cubic 
NLS from a 3D quantum many-body dynamic without a trap. The Elgart-Erdos-Schlein- 
Yau program consists of two principal parts: in one part, they consider the sequence of the 
marginal densities s 7^^'' [ associated with the Hamiltonian evolution e^'^^'^ipj^{0) where 



H 



N 



N 



-A. 



1 

N 



{ft -Tj)) 



l<:i<j^N 
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and prove that an appropriate limit of as — )■ oo solves the 3D Gross-Pitaevskii hierarchy 

k k 

(1.19) 2%(^) + 5^[A,„7(^)] =6o5^Tr,,^J<5(r,-r,+i),7^'+'^], for all A; > 1 . 



In another part, they show that hierarchy (1.19) has a unique solution which is therefore a 



completely factorized state. However, the uniqueness theory for hierarchy (1.19) is surpris- 



ingly delicate due to the fact that it is a system of infinitely many coupled equations over 
an unbounded number of variables. In ^H], by imposing a space-time bound on the limit 
of Ittv^I? Klainerman and Machedon gave another proof of the uniqueness in through 
a collapsing estimate originating from the ordinary multilinear Strichartz estimates in their 
null form paper [3H] and a board game argument inspired by the Feynman graph argument 



m 



Later, the method in Klainerman and Machedon [39] was taken up by Kirkpatrick, Schlein, 
and Stafiilani |37j, who derived the 2D cubic NLS from the 2D quantum many-body dynamic; 
by Chen and Pavlovic [9l [10], who considered the ID and 2D 3-body interaction problem 



and the general existence theory of hierarchy (1.19); and by X.C [16], who investigated the 



trapping problem in 2D and 3D. In [12] [13], Chen, Pavlovic and Tzirakis worked out the 



virial and Morawetz identities for hierarchy (1.19). In 2011, for the 3D case without traps, 
Chen and Pavlovic [Hj proved that, for (3 E (0, 1/4) , the limit of |77v'*| actually satisfies 
the space-time bound assumed by Klainerman and Machedon [39] as — )■ oo. This has 
been a well-known open problem in the field. In 2012, X.C. [TT] extended and simplified 
their method to study the 3D trapping problem for (3 G (0, 2/7]. 

The (3 = case has been studied by many authors as well [221 [3 UHl HSl [IB] • 
Away from the usage of the BBGKY hierarchy, there has been work by X.C, Grillakis, 
Machedon and Margetis [211 [321 [ISl [30] using the second order correction which can deal 
with e**^^'?/'^ directly. 

To our knowledge, this is the first direct rigorous treatment of the 3D to 2D dynamic 
problem. We now compare our theorem with the known work which derives nD cubic NLS 
from the nD quantum many-body dynamic. It is easy to tell that Theorem [1.2[ deals with 
a different limit than the known work [3l[2Il[23ll2ai25l[26l[2Zl[3Zl[l0l[l6l[IIl[I7] which 



derives nD NLS from nD dynamics. On the one hand. Theorem L2 deals with a 3D to 
2D effect. Such a phenomenon is described by the limit equation (1.16) and the coupling 
constant / |/ii(-2)| dz. The limit in Theorem 1.2 is with the scaling 



lim N^/uJ scat 



constant. 



instead of the scaling 



lim Nsca.t{N''^-^V{N^-)) = constant. 



in the known nD to nD work. 
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The main idea of the proof of Theorem 1.2 is to investigate the hmit of hierarchy (1.18) 



which at a glance is similar to the nD to nD work. However, in contrast with the nD to nD 



case, even the formal limit of hierarchy (1.18) is not known. 



Heuristically, according to the uncertainty principle, in 3D, as the 2;-component of the 
particles' position becomes more and more determined to be 0, the 2;-component of the 
momentum and thus the energy must blow up. Hence the energy of the system is dominated 
by its z-directional part which is in fact infinity as A^, — > oo. This renders the energy and 
thus the analysis of the x— component intractable. 

Technically, it is not clear whether the term 

o2 I 2 



00 



(k) 

tends to a limit as iV, a; — )■ oo. Since Tnuj ^ factorized state for t > 0, one cannot 



expect the commutator to be zero. Thus we formally have an oo — oo in hierarchy (1.18) as 



N,uj ^ oo. This is the main difficulty we need to circumvent in the proof of Theorem 1.2 



1.1. Acknowledgements. J.H. was supported in part by NSF grant DMS-0901582 and a 
Sloan Research Fellowship (BR-4919). X.C would like to express his thanks to M. Grillakis, 
M. Machedon, D. Margetis, W. Strauss, and N. Tzirakis for discussions related to this work, 
to T. Chen and N. Pavlovic for raising the 2D to ID question during the X.C.'s seminar talk 
in Austin, to K. Kirkpatrick for encouraging X.C. to work on this problem during X.C.'s 
visit to Urbana. We thank Christof Sparber for pointing out references |lj,|2j|. 

2. Outline of the proof of Theorem 11.21 

We begin by setting down some notation that will be used in the remainder of the paper. 
We will always assume a; > 1. Note that, as an operator, we have the positivity: 

-1 - dl +z^>0 

Define 

(2.1) 5/=^(l-A.^.+a;(-l-aJ + z|))V2 

We have Sj{(j){xj)h{zj)) = (1 — A^j0(xj) h{zj) and thus the diverging u parameter has no 
consequence when the operator is applied to a tensor product function (j){xj)h{zj) for which 
the 2;j-component rests in the ground state. 

Let Pq denote the orthogonal projection onto the ground state of —d"^ + z^ and P\ 
denote the orthogonal projection onto all higher energy modes, so / = Pq + -Pi? where 
/ : L^(]R^) — L^(]R^). Let Pq and P{ be the corresponding operators acting on L^(]R^^) in 
the Zj component, 1 < j < A. Then 

k 

(2.2) / = JJ(PJ+P/), where / : L2(m3A^) ^ ^^(M^JV) 

For a /c-tuple ex. = (cti, . . . , a^) with aj G {0, 1}, let Pq, = P^^ • ■ ■ P^'^. Adopt the notation 

\a.\ = ai + ■ ■ ■ + ak 
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This leads to the coercivity (operator lower bounds) given in Lemma A. 5 



We next introduce an appropriate topology on the density matrices as was previously done 
in lallEaEalElESlEniEZlEZllIOlini Denote the spaces of compact operators and 
trace class operators on Li^ (M^^) as ICk and respectively. Then (/Cfc)' = C\. By the fact 
that /Cfe is separable, we select a dense countable subset C Kk in the unit ball of 

/Cfc (so ||t/f'^''||op ^ 1 where |Hlop is the operator norm). For g C\, we then define a 

metric dk on C\ by 



4(7'", 7"=') 



1=1 



Tr J, 



(k) 



A uniformly bounded sequence 7^^''^ G C\ converges to 7'-*'^ G C\ with respect to the weak* 
topology if and only if 

hm 4(7!;,l,¥'=^) = 0. 

For fixed T > 0, let C ([0,T] , C\) be the space of functions of t G [0,T] with values in C\ 
which are continuous with respect to the metric dk- On C ([0, T] , C\) , we define the metric 

4(7^'U-),7^'U-))= sup 4(7^'U^),7('U^)), 

te[o,T] 

and denote by Tprod the topology on the space ®k^iC ([0,T] ,£j,) given by the product of 
topologies generated by the metrics dk on C ([0, T] , C\) . 

With the above topology on the space of marginal densities, we now outline the proof of 
Theorem 1.2 We divide the proof into five steps. 



Step I (Energy estimate). We transform, through Theorem 3.1, the energy condition (1.17) 
into an "easier to use" type energy bound in which the interaction V is not involved. Since 
the quantity on the left-hand side of energy condition (1.17) is conserved by the evolution, 
we deduce the a priori bounds on the scaled marginal densities 

k 

sup Tr n (1 - A.^, +uj{-l-d% + z|) ) 7!v'l ^ ^ 



supTrn(l-A.j7!J,U^' 



supTrP„7;j^^P^ < C'^w-il^l-il^l 



via Corollary 3.1 We remark that, in contrast to the nD to nD work, the quantity 

Tr(l-A,j7!;,l 

is not the one particle kinetic energy of the system; the one particle kinetic energy of the 
system is Tr (l — A^^ — cj9|J 7^^''^ and grows like uj. 

Step II (Co mpa ctness of BBGKY). We fix T > and work in the time-interval t G [0,T]. 
we establish the compactness of the sequence T]^^^^{t) = 



4.1 



In Theorem 

®k^iC ([0, T] , £^) with respect to the product topology Tprod even though there is an 00 — (X) 



N 



k=l 
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in hierarchy (1.18). Moreover, in Corollary 4.1, we prove that, to be compatible with the 
energy bound obtained in Step I, every limit point T{t) = {7'''^'*}^_-^ must take the form 

k 

7^^) {t, (x,, z,) ; (x'„ z'J) = ji'\t, X,; x'J J] h (z,) h {z'^) , 

i=i 

where 'j^^^ = Tr^^''-'^'^ is the x-component of 7'-*^^ 

Step III (Limit points of BBGKY satisfy GP). In Theorem 5.1, we prove that if r{t) = 
{^^'^''Ifeli is a ^ limit point of TN,uj{t) = W7v,L f with respect to the product 

~ 00 L ' J k=l 

topology Tprod, then {'j^^^ = Tr^ 7*-''''}^__^ is a solution to the coupled Gross-Pitaevskii (GP) 
hierarchy subject to initial data 7^^-* (0) = |0o) (^ol'^*' with coupling constant bo = J V (r) dr, 
which written in differential form, is 



'9tli'^ = E [-^^. ' ^i'T + ^0 E Tr.... Tr. [6 (r, - r,+0 , 7^'=+^)] 



Together with Corollary 4.1, we then deduce that {7^^^^ = Tr^ 7'-'^'*}^^ is a solution to the 
well-known 2D GP hierarchy subject to initial data 7^'^'' (0) = |0q) (^oI*^^ with coupling 
constant bo (J \hi {z)\ dz), which, written in differential form, is 

(2.3) za,7i^) = [-A.„7i'^] +&0 ( / \h {z)fdz) J2^T^... [Hxj-Xk+i),li'^'^] ■ 
j=i / 

Step IV (GP has a unique solution). When 7^ (0) = |0o) 

(00 1 ' know one solution to 
the 2D Gross-Pitaevskii hierarchy (2.3), namely 10) where (p solves equation (1.16). 

Since we have the a priori bound 



supTrn(l-A.,)7i'^^C^', 



the uniqueness theorem (Theorem 6.3) then gives that 'y^'^ 

N 



Thus the compact 



sequence r7v,aj(^) = |77vL| '^^^Y ^ ^ io'"(f^)+<^ limit point, namely 

' L ' J k=l 

k 

7^ = J]0(t,a;,)0(t,a;;.)/ii (z,) h,{z'^) . 
i=i 

By the definition of the topology, we know, as trace class operators 

k 

InI -^l[<P{t,x,)'^{t,x'^)h, (z,) h,{z]) weak*. 
i=i 

Remark 1. This is in fact the very first time that the Klainerman-Machedon theory applies 
to a 3D many-body system with (3 ^ 1/3. The previous best is (3 E (0,2/7] in [E] after 
the /3 G (0,1/4) work [TT]. Of course, we are not actually using any 3D Gross-Pitaevskii 
hierarchies here. 
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Step V (Weak convergence upgraded to strong). We use the argument in the bottom of p. 
296 of [27] to conclude that the weak* convergence obtained in Step IV is in fact strong. We 
include this argument for completeness. We test the sequence obtained in Step IV against 
the compact observable 

k 

i=i 

and notice the fact that ^T^vL) ^ 77vL since the initial data is normalized, we see that as 
Hilbert-Schmidt operators 

k 

^!v,l -^Yl<P{t,Xj)'(f){t,x'j)hi (zj) hi{zj) strongly. 
i=i 

Since Tr^yj^^^ = Tr7'^'^-', we deduce the strong convergence 



lim Tr 



via the Griimm's convergence theorem [501 Theorem 2.19] 



3. Energy estimate 

We find it more convenient to prove the energy estimate for ipjy^ and then convert it by 
scaling to an estimate for ipj^^^ (see (1.10)). Note that, as an operator, we have the positivity: 



-co - dl + w^z? > 



Define 



5^ = (l-A.,-u;-9j+c.^4 



2 2N1/2 



[I - UJ - Ar^ +^ 



2^2n1/2 



Theorem 3.1. Let the Hamiltonian he defined as in (1.7) with (3 G (0,2/5). Then for all 
£ > 0, there exists a constant C > 0, and for all uj,k ^ 0, there exists No^k, u) such that 

2 



(3.1) 



N,i. 



L2(M3iV) 



for all N ^ and all ^ e Ll {R^^) n V{H 



Proof. We adapt the proof of [211 Prop. 3.1] to accommodate the operator —u — dl. + u'^z'^ 
in place of —dl^ ■ The case /c = is trivial and the case k = 1 follows from the positivity of 
V and symmetry of We proceed by induction. Suppose that the result holds for k = n, 
and we will prove it for k = n + 2. By the induction hypothesis, 

{N-NU + H^^^r^^^) 
(3 2) " 
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For convenience, let 
Expand 

N / n 

l=n+l \l=l 

and substitute in both occurrences of the operator — Nu + Hn,u} in the right side of (3.2) 
to obtain four terms. We ignore the last (positive) one of these terms to obtain 

(3.3) (V^, {N-NUJ + i/^,^)"+V) > C"iV"(I + II + III) 
We have 

N n 
il,h=n+l j=l 

In this double sum, there are (A^ — n){N — n — 1) terms where ii ^ £2 that are all the same 
by symmetry, and there are (A^ — n) terms where ii = £2 that are all the same by symmetry. 
We have 

n+2 n+1 

(3.4) l = iN-n){N-n-l){^,l[ S]^) + {N - n){^, J] S]^) 



the first of which will ultimately fulfill the induction claim. In (3.3), we also have 

N n n N n 

II + III = 2 j2 E n ^pi"^) + E n 

^i=n+1^2=l i=l e=n+l j=l 

N n 

e=n+l j=l 

Exploiting symmetry this becomes 

n+l n+1 

(3.5) II + III = 2(Ar - n)n{ij, J] S^tP) + 2{N - n) Re{ij, J] S^H'^^^ij) 

j=i j=i 

In the first term, we have applied the permutation that swaps ii and n+1 and £2 and 1. 
In the second and third terms, we have applied the permutation a that swaps i and n + 1. 
Strictly speaking, this permutation maps Hj^^ to Hj^^^ where 

l<i<j<N 

where ±1 is chosen according to the affect of the permutation on the pair The dis- 

tinction between Hj^^ and H}^ ^^^^ is inconsequential for the remainder of the analysis (and 
in fact Hj^^ = Hj^^^^^ if V is even), so we have ignored it in (3.5). The first of the terms 
in (3.5) is positive - it is the second term that requires attention; in particular, we have to 
manage commutators. 
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Assuming > 2n + 2, we substitute (3.4), (3.5) into (3.3) to obtain 

n+2 



n+1 



(3.6) 



i=i j=i 



n+l 



-2C"iV"(iV - n) Re{^, JJ S]Hj,^^^) =: D + E + F 



The first two terms, D and E, in (3.6) are positive. The third term F will be decomposed 
into components, some of which are positive and others that can be bounded in terms of the 
first two terms appearing in (3.6). In the expression for Hj^^, there are 

• ^(n + l)n terms of the form V{ri — Vj) for l<z<j<'n,-|-l. 

• (n + 1) (A^ — n—1) terms of the form V{ri — rj) ioi 1 < i < n + 1 and n + 2 < j < N. 

• ^{N — n — V){N — n — 2) terms of the form V{ri — r^) for n + 2 < i < j < A^. 

For convenience, let 
Using symmetry, we obtain 

n+l 

F = 2C''N"\N -n){n + l)nRe{^, JJ S]Vi2ij) 

n+l 

+ 2C"Ar"(Ar -n){n + l){N -n-1) Re{ij, J] S]V^n+2)i^) 



i=i 

n+l 

+ C'-'N^iN -n){N-n-l){N-n-2) Re{ij, J] 5|V(„+2)(n+3)V') 

i=i 

=: Fi + F2 + F3 

The last term F3 is positive since each 5*^ for 1 < j < n + 1 commutes with V(„+2)(n+3)- We 
will show Fi > —\E and F2 > —\D provided A^ > No{n), which together with (3.6) will 
complete the induction argument. We have 

n+l 



Fi = 2C"Ar"(Ar - n)(n + l)nRe(^, JJ S]Vi2^P) 



i=i 



2C"Ar"(Ar - n)(n + l)nRe / (/, SfSiVi2f)rur, dr^--- drN 

Jr3,...,rM V ' 



■■■■Fi 



where / = YYS=3 Sji^- We can regard r^, . . . ^r^ as frozen in the following computation, so 
to prove < ^E, it will suffice to show that 



(3.7) 



F,| < Kll Will, 
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Toward this end, we have 

\h\ = \{Slf, VuSif) + 2(5?/, v., 1^12 ■ VrJ) + {S'J, {AM f)\ 



< 



hfh'i^Le:,^ ||^12||l-L3^ II 'S's / || L2^L2^ 



+ \\Slf\\^^L<^J\Vr,V, 



12 



Too r 3/2 



IV.JII 



By evaluation of 



1211 roo r6/5 



J 2 r oo 



II ^12 II L3^ ~ 



2/3-1 



l|V.,K 



12|| r3/2 



2/3-1 



12\\ r6/5 



/3-1 



the above estimate reduces to 

\F,\ < (iVa;V^)^^-iS,V||LP^LBj|52Vl|L2^L2^ + (iV^^/')'^-l^?/l|L2^Lo;J|V.J|U2^^a 



+ (iVu;V2)t^-||5?/|U2,a 



L2 



Applying Lemma A. 4, this reduces further to 

\h\ < (Arc.V^)^^-^c.V6||5,2^2/IU.^L2 

+ (Aru;V^)^^-u;VV/iS?S2/|U2^.2j|52VlU,.2 

+ (Aru;V2)|/3-l^l/6^1/4||^2^^^|| . . „c2 



1^2 <> ii-^ri-^r2 



Hence we need /3 < ^ and conditions (3.13), (3.11) below to achieve (3.7). 



Let us now establish F2 > —\D. We have 



n.+l 



F2 = 2C"Ar"(Ar - n)(n + l)(Ar - n - 1) Re(^/., J] '5,Vi(n+2)^) 

i=i 



2C"iV"(iV - n)(n + l)(iV - n - 1) / (/, ^iVi(„+2)/)n,r„+2 c?^2 ■ ■ ■ c^r„+ic^r„+3 ■ ■ ■ 



■■Fi 



where / = IlS ^j'^- Now 



r'n + 2 



~ '^(Z, ^l(n+2)/)r-ir„+2 r'ir„+2 

+ (C^zi/j V'i(„+2)C^zi/)rir„+2 + (/) ^^^l/)rir„+2 
= : -^2,1 + -^2,2 + -^2,3 + -^2,4 

Note that ^2,3 and F2,4 are positive and can thus be disregarded. To prove F2 > —\F), it 
suffices to prove 

(3.8) 

But 



|F2,l| + |F2,2|<^n-l||5iS„+2/||i2 ,2 



•■l ''n+2 



By Lemma 
(3.9) 



A.4 



|^2,l| ^t^ll/l|L2 ||\4(„+2)||^^^3/2 

and ||\/i(„+2)Loor3/2 ~ (iVa;^/2)/3-i^ obtain 



■^^l-^^n + 2 
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The upper bound in ( 3.8[ ) will be achieved provided (3.12) below holds. Also, 

1-^2,21 
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^21 Vl(n+2)|| roo r3/2 



< 1 1 / 1 1 L2 LP 1 1 ^21 "^4 (n+2) 1 1 ^3/; 

2/3-1. By Lemma 



A.4 



Note that ||(92jVi(„ 

and ||/||l2 re < w^/^||S'„+2/||l2 L2 , • From this, it follows that 



'■l '■n + 2 



(3.10) 



2,2 



The upper bound in (3.8) will be achieved provided (3.13) holds. By (3.9), (3.10), we obtain 
(3.8), completing the proof. Let us collect the conditions on and u. We have 



5 fl_ 1 



(3.11) 


(Aru;i/2)l/^-V/i2<<n-2 




(3.12) 




3/3+1 1 

<^=^ A^ > w 1-/3 ni-/5 


(3.13) 


(iVu;i/2)2/3-i^5/6 ^ ^-1 





The requirement that (3.11), (3.12), and (3.13) hold is imposed in the definition (1.15) of 

v{(3). □ 



Now consider the rescaled operator (2.1) so that 



(S'jV')(t,XAr,ZAr) = UJ^ /^{S jll)){t,yiN , VuJZn) ■ 



3.1 



into statements about ^, Sj, and 7^'*^ 



We will convert the conclusions of Theorem 
that we will then apply in the remainder of the paper. 

Corollary 3.1. Let i'N^^it) = e^^-^^'^-^^YwlO) '^''^d {TtvLI'^)} marginal densities as- 

sociated with it, then for all to > 1 , k > 0, N > co'"^^^'^'^, we have the uniform-in-time 
bound 



i=i 



(3.14) 

Consequently, 

(3.15) Trn(l-A„j7!;.! 



L2 



n(i-A.,)i/^^^,jt) 



L2 



and 
(3.16) 



Af,wllL2 



Proof. Substituting (1.10) into (3.1) of Theorem 3.1 and rescaling, we obtain 

2 



(3.17) 



n ^j^N,u 



L2 



18 XUWEN CHEN AND JUSTIN HOLMER 

Since N — Hn,ui — Nu is self-adjoint and [Hn,u), N — Hn,u) — Nu] = 0, 



Hence by (|3.17|) 



L2 



where the last estimate follows from the hypothesis (1.17) of Theorem 1.2 



The inequality (3.15) follows from (3.14) and 



follows from (A. 29) and (3.14). By Lemma 



inequality on the right of (3.16) follows by Cauchy-Schwarz. 



A.6 



A.27). The inequality on the left of (3.16) 

TrP„7!JL^/3 = {Pcx^N,u,Pf3^NJ, so the 



□ 



4. Compactness of the BBGKY sequence 
Theorem 4.1. The sequence 



which satisfies the cxd — oo BBGKY hierarchy (1.18), is compact with respect to the product 

topology Tprod- For any limit point T{t) = {7'''^^}^_-,^ , 'j^''^ is a symmetric nonnegative trace 
class operator with trace bounded by 1. 



We establish Theorem 4.1 at the end of this section. With Theorem 4.1, we can start 



talking about the limit points of rjv^^(t) = {T^l^j^Li 



Corollary 4.1. Let r{t) = {7^'^-'}^;^ be a limit point ofrN,Lu{t) = {7iv,L}fc^i '"'^^^ respect to 
the product topology Tprod, then 7'-'^^ satisfies 



(4.1) 



TrJ](l- A,j7('=) ^C'^ 



(4.2) 



7^'=) (t, (x,, z,) ; (x'„ z'J) = t'\t, X,; x'J J] h {z,) h, {z'^) 

i=i 



Proof. The estimate (4.1) is a direct consequence of (3.15) in Corollary 3.1 and Theorem 



4.1 The formula (4.2) is equivalent to the statement that if either o; 7^ or /3 7^ 0, then 



Pccl^^'^Pp = 0. This is equivalent to the statement that for any J^''^ e ICk, Tr J^^^ P^j^^'^ Pf3 
0. However, 



(4.3) 

By Lemma [A.6 



TrjWp«7(^)p^= hm Tr j(^)p«7^),P^ 

(A'',a;)->oo ' 
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and by Cauchy-Schwarz and (3.16) 



Tr jWp„7!J;li^/3| < ||J('=)||op||P«^^,J|L2(M3^)||P/3t^^,J|L2(R3.v) <C'=a;-i 



5l"|-|l/3| 



Hence the right side of (4.3) is 0. 



□ 



compactness of 7jv|^ for fixed k with respect to the metric dk- By the Arzela-Ascoh theorem, 
this is equivalent to the equicontinuity of 7!^''^, and by [271 Lemma 6.2], this is equivalent 
to the statement that for every observable J^'^^ from a dense subset of }C{L'^ (R.^^)) and for 
every e > 0, there exists 5{J^''\e) such that for all ^1,^2 £ [0,7"] with \ti — ^2! ^ S, we have 



Proof of Theorem 4-1 ■ By the standard diagonalization argument, it suffices to show the 



(4.4) 



sup 

N,u] 



We assume that our compact operators J^^^ have been cutoff as in Lemma A. 7 Assume 



ti < t2. Inserting the decomposition (2.2) on the left and right side of Tnuj^ we obtain 

~{k) 



where the sum is taken over all fc-tuples a and f3 of the type described above. 



(4.5) 



To establish (|4.4|) it suffices to establish, for each ct and f3 



sup 



Below, we establish the estimate 



Tr jWp«7!J,l^/3(t2) - Tr /'=)P«7!f,l^/3(ti^ 



(4.6) 



< It, 



if both CK = and /3 = 



lmax(l,(jj^ 2!°! 2 1^1) otherwise 



Estimate (4.6) suffices to prove (4.5) except when \ct\ = and |/3| = 1 or vice versa, in which 
case it yields the upper bound u^^'^\t2 — ti| with the adverse factor u^^"^. On the other hand, 
we can also prove the (comparatively simpler) bound 



(4.7) 



Trj('=)p„7!vV/^(t2)-TrjWp^7!vVM^i)l^^"^'"''^'''' 



that provides no gain as ^2 ~^ ti, but a better power of 00. By averaging (4.6) and (4.7) in 
the case |q;| = and |/3| = 1 (or vice versa), we obtain 



which suffices to establish (4.5) 



20 



XUWEN CHEN AND JUSTIN HOLMER 



Thus, it remains to prove both (4.6) and (4.7), and we begin with (4.6). Hierarchy (1.18) 
yields 

k k 

(4.8) idtP^^f^Pp = \-A,^, Po7S{,l^/3l + + 'I Po.l%lPp 



Let 



1 



+ 



N-k 



N 



Ph 



I = -zETrjW[-A.^,P„7SvV/3] 



(4.9) 



II = -u^ Y Tr J''^ l-dl + zl P^if^M 
III = -ziV-i Y Tr /'•)Pc.[V5v,.(ri-r,),7SJ,l]i'/3 

l<i<i<fc 
k 



IV = -i'- 



N-k 



N 



Then it follows from (|4.8|) that 
(4.10) 



dt Tr J('=)p«75f>^ = I + II + III + IV 



First, consider I. Applying Lemma A. 6 and then integration by parts, we obtain 

k 

I = ((jWA,^.P„^,P^V^)r, - {J^'^Paij,P^A,^^),^) 

= ((j(^-)A.^P„V^,P;3^)r, - (A.,.jWP«V^,P^^)rJ 

i=i 



Hence 
(4.11) 



1 1 op 



|A^^.J^'')||op)||PaV^||L2(K3iV)||P^7/'||i2(K3iV) < Cfc_j(fc) 



|I|<E(||J«A, 
i=i 

where in the last step we applied the energy estimate. 
Now, consider II. When ck = and /3 = 0, we use that 

k 

II = -uz Y Tr J^'^ [1 - 9J + zl Pcl^Pf^] = 
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Otherwise, we proceed directly from (|4.9|), applying Lemma A. 6 and integration by parts to 

k 



obtain {Hj = —d"^. + 



k 



cut 



Hence 



By the energy estimates, 



(4.12) 



Now, consider III. 



II 







if a = and /3 = 



< Cj.j(k) uj^ al"! 2 1^1 otherwise 



l<i<j<k 

= -tN-' J2 {j''''^PcyNM-r,)i^,P(ii^) - {Pc.^,J^''^P(iVM,^{r,-rj)^) 

l<i<j<k 

Let Li = (1 - 



Then 



III = -iN-' {J^^'^Pc.ULjWijULji,, P^^) - (P„^, J^'^^PpUL.WijUL,^) 

l<i<j<k 



Hence 



|in| < iV-i||jWL,L,||opl|lV,,||op||L,L,V^|U2(R3iV)||P;3Vl|L2 

+ N-^\\P^^\\L2^^3N)\\J^^^LiLj\\,jW,j\\op\\LiLj^\\L2 



By Lemma A.l, ||VFjj||op ^ ||VAf,a;||Li = (independent of N, uj), and hence the energy 

estimates imply that 



(4.13) 

Now consider IV. 



IV = -i- 



N-k 
N 



|III| < 



J2 {{J^^^PcVNArj - n+i)iJ,P^ij) - {j'^'^Pc.i^,P,3VNAr, - rk+iW) 
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Then, since J^^^^L^+i = Lk+iJ^^\ 



IV = -i 



N -k 
N 

N-k 



k 



— l- 



N 



^{LjJ^^^P^Lk+ii),Pi3Wj^k+i)LjLk+i^) 



Estimating yields 

k 

|IV| < $^(||j('^^,||op+ ||i.,/'^)||op)||W^,(fc+l)l|op|| LjLk+li'\\L'^(K:-^r^) \\Lk+l'4'\\L'^ 

i=i 



By p.lSD , 
(4.14) 



IIVI < C, 



k,jik) 



Integrating (4.10) from ti to t2 and applying the bounds obtained in (4.11), (4.12), (4.13), 
and (14. 141), we obtain d^el). 



Finally, we proceed to prove (4.7). We have, by Lemma [A. 1 



< 2 sup I (JC^^P,^^, Jt), P^iJj,Jt))r, I 
t 



where in the last step we applied (|3.16|). 
According to Corollary 4.1 



□ 



the study of the limit point of Tjy^^^it) = {t^I^I is directly 

N ' I ' J fc=l 

related to the sequence r,,jv,c.(t) = (tI^^L = Tr. tSJI) e (Bk^iC ([0, T] , £i (M^'c)) . We 

L ' ' 'J k=l 

will do so in Section |5j We end this section on compactness by proving that Tx,N,co{t) is 
compact with respect to the two dimensional version of the product topology Tprod used in 
Theorem 4J^ This proof is not as delicate as the proof of Theorem AA_ because we do not 
need to deal with oo — oo here. 



Theorem 4.2. The sequence 



it) 



~{k) 
'~1x,N,ui 



is compact with respect to the two dimensional version of the product topology Tprod used in 
Theorem\4-l ■ 
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I I (k) 

Proof. Similar to Theorem |4.1[ we show that for every observable Jx from a dense subset 
of /C (L2 (M2fc)) and for every e > 0, 36{ji''\e) s.t. Vti, ts e [0, T] with \ti - tsl ^ 5, we have 



sup 

N,uj 



^ e. 



We utilize the observables Jx'elC (L^ (M^*')) which satisfy 



< oo. 



op 



4.1 



~ik) 
since -fl ' 



Here we choose similar but different observables from the proof of Theorem 
acts on (M^'^) instead of (M^'^) . This seems to make a difference when we deal with 
the terms involving 7^|^ or 7'-'^''. But ji'^'' does nothing on the z variable, hence 



lop 



((v.,> + 5.,) 



op 



((v.,> + s.,) 



op 



op 



lop 



I.e. 



\LjJx ""^llop; 11-^1 ^ Jx ''^jllop; \\LiLjJx ^ ^Lj ^||op sud \\L^ ^Lj ^ ji ^ LiLjWqp are all fi- 



nite. It is true that ji^^ and the related operators listed are only in (L^ (M'^'^)), but this 
is good enough for our purpose here. 

Taking Tr^ on both sides of hierarchy (1.18), we have that "J^^l^^ satisfies the coupled 
BBGKY hierarchy: 

(4.15) = X:[-A.,.,CU+4ETr.K.(r.-r,),7!;i 



i<j 



+ 



A^- A; 



Tr^,+i Tr^ Vn,u^ (rj - r^+i) , 7$^ 



~(fc+i) 



Assume ti ^ ^2, the above hierarchy yields 



Tr J?) 



-A ^^''^ 



i<j 



t2 



Tr 



A^- A; ^ 



AT 



t2 



Tr J?) 



t2 



i<j 



N-k 



r,(fc+l) 



AT 



t2 



i=l ^1 
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For I, we have 



Tr Jl'^) 



-A ^^''^ 



(Vxj)^ yl'xNui (1 commutes with everything) 



Tr (V.^,)-^ (V.,y 7?i,. (v.,) - Tr (V.^.> {V.y (V.,> tI^. (V.,> 
(p.,r4'^{V.,)\\ +||(V.,>4'^)(V.^r1| )Tr(V.^,>7ltL(V.^.> 

V op op/ 



^ C,Tr(V.,y7SJ,l 
^ Cj (Corollary [3l|. 

for II and III, we have 



II 



Tr Ji.'^) 



J 1 1 op 



and similarly, 



III 



Tr J?) 



Tt LJ^ L'^l-^^J^''^ LjLk+iWj(^k+i)LjLk+i^^^^^^ LjLk+i 

J 1 1 op I 



^ (||LtV(^)L,||^p + ||L,Ji^)LTi|| J ||iy,(.+i)|LpTr L,W7!vU'^L, W 



Up to this point, we have proven uniform in time bounds for I - III, thus we conclude the 



compactness of the sequence Tx,N,u]{t) = |7i^Ar,aj| 



□ 



5. Limit points satisfy GP hierarchy 

Theorem 5.1. Let T{t) = {7^^'^}^^i be a N ^ Umit point ofTN^^{t) = {tSJ.I}^ ^ 

with respect to the product topology Tprod, then {7^'^^ = Tr^ 7*-'^'*}^-^ is a solution to the coupled 
Gross-Pitaevskii hierarchy subject to initial data 7^'^-' (0) = |0o) (0ol'^'^ '"^^^^ coupling constant 
bo = J V (r) dr, which, written in integral form, is 

(5.1) 7i'^ = U^'\t)^i'^ (0) - z6o f ^^'^^ - Tr.,,, Tr, [5 (r, - r,+0 , ^^'^^^'^ (.)] ds, 

3=1 
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where 



We prove Theorem 5.1 below. Combining Corollary 4.1 and Theorem 5.1, we see that 
7^'^'' in fact solves the 2D Gross-Pitaevskii hierarchy with the desired coupling constant 
bo {J\h {z)\Uz). 

Corollary 5.1. Let T{t) = be a N ^ ^u^C^+e Umit point o/rjv,c.(t) = {inIV 

^ I. ' J k=l 

with respect to the product topology Tprod, then {7^^-' = Tr^ 7''*''^}^__|^ is a solution to the 2D 
Gross-Pitaevskii hierarchy subject to initial data ^^J'^ (0) = |0q) (0oI'^'^ '"^^^^ coupling constant 
bo (/ 1^1 dz) , which, written in integral form, is 
(5.2) 

= f/W(t)^W (o)_,6„ (^j \h, [z)t dz^ [ U^'\i-^) Tr..,. [5 (x, - xu+i) , li'^'^ (s)] ds. 

Proof. We compute the k = 1 case explicitly here. Written in kernels, the inhomogeneous 
term in hierarchy (5.1) is 



ibo J U^^\t-s)ds J 6{zi - z[)dzidz[ J S{ri - r2)^'^'^\ri,r2,r[,r2)dr2 
-ibo / U^^\t-s)ds / 6{zi - z[)dzidz[ / 6{r[ - r2)j'^^\ri,r2,r[,r2)dr2 



which, by Corollary 4.1 



IS 



= ibQ J U''^\t - s)ds J S{zi - z[)S{ri - r2)'yf\xi,X2,x[,X2) 

xhi{zi)hi{z2)hi{z[)hi{z2)dr2dzidz[ 
-ibo J U^^\t - s)ds J 5{zi - z[)5{r[ - r2)^^^\xi,X2,x[,X2) 

xhi{zi)hi{z2)hi{z[)hi{z2)dr2dzidz[ 

Further simplifications lead to 

= ibo Ju^^\t~s)ds J 6{xi - X2)^f\xi,X2,x[,X2)\hi{zi)\'^dx2dzi 

-ibo / U^^\t - s)ds / 5{x[ - X2)^f\xi,X2,x[,X2)\hi{z[)\'^dx2dz[. 



In summary, we have 

tbo J U^'\t-s)TT^,TT, [6 (ri-r2),7^2) (s)] ds 

= ib 



'0 (^j \h, {z)tdz^ j f/(2)(t-s)Tr,, [5(xi-a;2),7l'^ {s)] ds. 



□ 
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Proof of Theorem 5J_. By Theorems |4.1[ |4.2[ passing to subsequences if necessary, we have 
hm supTrj('=)f7;;)^(t)-7«(t)) = 0, V Z'^) G ^ (L^ (M^fe)) 

(5.3) 



AfjO;— j-oo f 



hm supTrjW(7S^,,(t)-7i'^)(t))= 0, V G /C (L^ (M^'^)) . 



We estabhsh (5.1) by testing the hmit point against the observables Jx E K, (L^ (M^'^)) as 
in the proof of Theorem 4^ We wiU prove that the hmit point satisfies 

(5.4) TrjW7i'Uo) = Tr4'=)|0o)(0or' 

and 

TrJi'^Hi^) (t)=TrjW^(^) (t)^^^^ (0) 

(5.5) 



-ihoi^l Tr f/('=)(t - s) [5 (r, - r^+i) , 7^'=+^) {s)] ds. 
,=1 -^0 



To this end, we use the coupled BBGKY hierarchy (4.15) satisfied by 7^^^]^^, which, written 
in the form needed here, is 

Tr4^)7SL(t)=Tr4^)f/(^)(t)7it^,,(0) 

- ^ E 4'^^^'^ it - s) \V^^^ (r, - r,) , ^f^^ {s) ds 



— I 



i<j 



N-k 



N 



J2 Tr (t - s) [V5v,. (r, - r,+i) , 7!?:'^ (^) 



By (5.3), we know 



hm TrjW7i':^,^(t) = Tr 4'=)7i'=) (^) , 
hm TrjWf/W(t)7l'^_^(0) = Tr (t) 7!^) (0) . 

By the argument that appears between Theorem 1 and Corollary 1 in [42j, we know that 
assumption (b) in Theorem |1.1[ 

^ivL (0) ~^ 100 ® ^1) (00 ® hi\ , strongly in trace norm , 



in fact implies 



77vL (0) ~^ 100 ® ^1) (00 ® hi\^^ , strongly in trace norm . 
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Thus we have tested relation (5.4), the left-hand side of (5.5), and the first term on the 



right-hand side of (5.5) for the limit point. We are left to prove that 



lim — 

N,aj-^oo JM 



lim [1--]D 



bo f Tr 4'=)f/«(t - s) [5 (r, - r^+i) , 7^'+'^ {s)] ds. 
Jo 



First of all, we can use an argument similar to the estimate of III and IV in the proof of 



Theorem 4.1 to show the boundedness of \B\ and \D\ for every finite time t. In fact, noticing 



that U^''^ commutes with Fourier multipliers, we have 







\B\ ^ 


f 




Jo 



Tr {t - s) [V^,. (r, - r,) , ^f^, {s)\ ds 

j ds\Ti Lt'Lj'J^'^L.L.U^'^ [t - s) W.jLiLj^^^l (s) L,L, 
- Tr L,L,4^)LriLTif/W {t - s) UL.^f^^ {s) L,L^W,,\ 



^ Cjt. 



Hence 



hm — 

N,uj->-oo JM 



kD 

lim — — = 0. 

Af,cj-s>oo N 



To prove 
(5.6) 



lim D 



Tr 4'^U^'\t - s) [6 (r, - r,+i) , 7^'+') (s)] ds, 



we need Lemma A. 2 (stated and proved in Appendix |A]) which compares the 5— function 
and its approximation. We choose a probability measure p E (M^) and define p„ (r) 



(^) .In fact, p can be the square of any 3D Hermite function. Write Jg% = Jx'^^U^^^ {t — s) 
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we then have 

Tr it - s) (r^,. (r, - r,+0 7!?;'^ (^) - boS (r, - r,+0 7^'=+^) (s)) 

Tr Ji^^i (VW,a; [rj - Tk+i) - bo6 (r^ - rt+i)) 7!f,t^^ (s) 
Tr jS (5 (r^ - Tfc+i) - p„ (r^- - r^+i)) 75{+^^ (s) 
Tr JSp„ (r, - r.^i) (7!;:'^ (^) - 7^'+^^ (^)) 



I + II + III + IV 



j__ 1 

We take care of I first because it is a term which requires jy_> u^p 2, Write 
:^V{x, we have Vn^^ = {N^/uJ)^^ V^{{N^/^u)^ r), Lemma A. 2 then yields 



I ^ 



Cbo 



{N 



x(ll^.4%iop+ll^7'4%IL,)TrL,L 
^ ^ (/K.(r)|rrdr) 



■fc+i77v,t^^ (■5) LjLk+i 



Notice that (/ VL(r) |r|'^(ir) grows like (y/u)'^, so / ^ Cj I ^ — j which converges to 

_L 1 I 

zero as A^, w — )■ 00 in the way that N ^ uj'^f^ . More precisely, 



lim / = 0. 



So we have handled I. 
For II and IV, we have 

II ^ C'6o«'^(||L,4'=)LTi||^^+||LTV(^U,J|jTrL,W^^^^^^ (Lemmag^ 

^ Cja"" (Corollary [3l| 
IV ^ Cb,a-(^\L,J^^L-:\^+\\L-'j^^L,\l^TTL,^^^^^ (L 
^ Cja"" (Corollary [4l| 

which converges to as a — )■ 0, uniformly in A^, 
For III, 



emma A. 2) 



Ill < &0 

+&0 



Tr JSp„ (r, - r,+0 —±— i^^^l'^ (s) - (.)) 
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The first term in the above estimate goes to zero as A^, w — i- oo for every e > 0, since we have 
assumed condition (5.3) and J^'l^tPa ('^j ~ '''fc+i) (1 + ^^-^fc+i)^^ is a compact operator. Due to 



the energy bounds on 7^^^"'^'' and the second term tends to zero as e — 0, uniformly 

in A^. 



Combining the estimates for I-IV, we have justified hmit (5.6). Hence, we have obtained 
Theorem 15. 1[ 

□ 



6. Uniqueness of the 2D GP hierarchy 

For completeness, we discuss the uniqueness theory of the 2D Gross-Pitaevskii hierarchy. 
To be specific, we have the following theorem. 

Theorem 6.1 ([16, Theorem 3]). Define the collision operator Bj^k+i by 

Suppose that \ 7!^'* \ solves the 2D constant coefficient Gross-Pitaevskii hierarchy 
I J k=i 



(6.1) 



subject to zero initial data and the space-time bound 



(6.2) 







k 

n 



dt ^ 



for some C > and all 1 ^ j ^ k. Then V/c, t G [0, T], 

k 



n 



0. 



L2 



Proof. This is the constant coefficient version of flGl Theorem 3]. W. Beckner obtained 
the key estimate of this theorem independently in [5]. Some other estimates of this type 
can be found in [lH [29]. K. Kirpatrick, G. Staffilani and B. Schlein are the first to obtain 
uniqueness theorems for 2D Gross-Pitaevskii hierarchies. One will find their Theorem 7.1 in 



by replacing |V|2 by (V)2''''^ in the statement of the above theorem. 



□ 



To apply Theorem |6.1| to our problem here, it is necessary to prove that both the known 
solution to the 2D Gross-Pitaevskii hierarchy (namely |0) , where solves the 2D cubic 
NLS) and the limit obtained from the coupled BBGKY hierarchy (4.15), satisfy the space- 



time bound (6.2). It is easy to see that 



verifies the space-time bound (6.2) because 



it is part of the standard procedure of proving well-posedness of the 2D cubic NLS. We use 



the following trace theorem to prove the space-time bound (6.2) for the limit 
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Theorem 6.2 ([3T, Theorem 5.2]). For every a < 1, there is a Ca > such that 



(fc+l) rl ( r2 (Ta,2k\ 



^C.Tr n(l-A.,) 



for all nonnegative 7! ^ ^ G £^ (L^ (M^'')). 

We can combine the above theorems so that it is easy to see how they apply to our problem. 
Theorem 6.3. There is at most one nonnegative operator sequence 



that solves the 2D Gross- Pitaevskii hierarchy (6.1) subject to the energy condition 

7. Conclusion 

In this paper, by proving the limit of a BBGKY hierarchy whose limit is not even for- 
mally known since it contains (00 — 00) , we have rigorously derived the 2D cubic nonlinear 
Schrodinger equation from a 3D quantum many-body dynamic and we have accurately de- 
scribed the 3D to 2D phenomenon by establishing the exact emergence of the coupling 
constant (^J \ hi(z)\ dz). This is the first direct rigorous treatment of the 3D to 2D dynamic 
problem in the literature. 

Appendix A. Basic operator facts and Sobolev-type lemmas 
Lemma A.l (1211 Lemma A.3]). Let Lj = (l - A,.J ^ ^ y/^g^ /^^^g 

IlLr^LT V (r, - r,) Lr^Lj\^ ^ C \\V\\^, . 

Lemma A. 2. Let p E (R.^) be a probability measure such that J^^ (r)^ p (r) dr < 00 and 
let (r) = a^^p (^) . Then, for every k G (0, 1/2) , there exists C > s.t. 



|Tr J(^) (p, (r, - rfc+i) - 5 (r, - r^+O) 7^'+'^ | 

^C(^j p{r) \r\-dr^ (ll^.-^^'^^/'ILp + 1 1 ^7 1 1 op) Tr L.-L.+iT^'"''^^. W 



for all nonnegative 7('=+i) g C} (L^ (]R3''+3)) . 

Proof. We give a proof by modifying the proof of [371 Lemma A. 2]. We remark that the 
range of k, is smaller here because we are working in 3D. It suffices to prove the estimate 



for k — 
write 
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1. We represent 7^^) by 7^^) = ^. Xj \(p.) ((^^.|, where (pj e (R^) and ^ 0. We 



3 

= (^J' iPo^ i^^ -ri)-8 (n - rs)) ifj) 



where ip^ = (J*^^^ ® l) (fj- By Parseval, we find 



lii'j, - ^2) - 5(ri - r2))<^j)\ 

^ / \i^ji^i,QmCi,Qmi + ^2-^'i-Q\ y"p(r)(e--(^-«i)-l)dr|deicie2« 



Using the inequahty that G (0, 1) 



^ «'^irriei-e'ir 



we get 



',->(Pa(n -r2) -(5(ri -r2))(^.)| 



p(r)lrrdr) j e2)ll^,(e;, ^2)1^(^1 + C2 - - r2)c^^W^2C?CW^^^ 

+ a^{j p{r)\r\^dr) j e2)ll<^,-(e;, e2)l<^(ei + ^2 - - e;)^^ei^^e2<irfe; 



= q;'^( y p(r)|r|'^cir)(I + II). 
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The estimate for I and II are similar, so we only deal with I explicitly. 



i'ii^i,^2) d^^d^2d^'id^2 



+ - / 5(^1 + ^2-^1-^2) — 



d^id^2 



1 



1 

e 

^ e 



1 



(e'l + e'2 - 



, 2(l-«;) 



{ipj, LlLlipj) sup / \27\2dv + 7 {^j, LlLl^j) sup / 



1 



When K e [0,1/2), 



sup 



1 



1 



dr] < 00, 



sup 



2 



and hence we have (with e = \\LiJ^^^Li^\\-^), 
|IVJ« (pjri-r2)-5(ri-r2))7('+')| 

^ C p(r) |rr dr^ a'^ Tr J^L^L^ 7(^)7(2) + 1 Tr L2l27(2)^ 

<C(^I p{r) \rrdr^ H^^^^^^^^i ILp \\LiJ^'^L^\, + ^) TrL?L^7(^) 

□ 

Lemma A. 3 (some standard operator inequalities). 

(1) Suppose that A>0, Pj = P* , and I = Po + Pi. Then A < 2PoAPo + 2PiAPi. 

(2) If A>B> 0, and AB = BA, then A" > 5" for any a > 0. 

(3) // Ai > A2 > 0, Bi > B2 > and AiBj = BjAi for all 1 < i,j < 2, then 
AiBi > A2B2. 

(4) IfA>0 and AB = BA, then A^/^B = BA^^. 

Proof For (1), pi/^/f = \\A^^^{Po + Pi)f\\^ < 2\\A^/^Poff + 2\\A^/^Piff. The rest are 
standard facts in operator theory. □ 
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Recall that 

Lemma A. 4 (Estimates with cj-loss). Suppose f = f{x,z). Then 
(A.l) l|V./|U.<a;V2||5;||^, 
(A.2) \\fU<^'^'\\Sfh. 
(A.3) ||V./|U6<a;2/i5VlUp 
(A.4) \\f\\L?^<u^'^'\\S'f\\^ 

The factors of u appearing here are seen to be optimal by taking f{x,z) = g{x)h^{z), 
where g{x) is a smooth bump function. Then S"^ f = (1 — /S.x)g{x)h^^{z) and hence 

\\Sf\\% = (S'fJ) = {{1-A^)gh^,h^) = {\\g\\% + ||V.^?||i.)||/.J|i. 

which is ^-independent. Also, 115*^/11^2 = ||(1 — Ax)g\\Li independent. On the other 
hand, it is apparent that \\Vrf\\^ = uj^^^, = W^rfh^^ = ^'^'^ and = 

which demonstrates sharpness of the estimates. 

Proof. Recall J = Pq + -Pi- First, we establish 

(A.5) llV.Pi/IU^ < II5/IU2 

(A.6) ||Pi/|Ua < 

(A.7) ||V.Pi/|U6 < II^VIU^ 

(A.8) ||Pi/|U^ < II^VIU? 

Note that PjS"^ = S'^Pj. By the definition of S, 

Pi(l - A, + uh^)Pi = S^Pi + uPi 
By spectral considerations 2ujPi < PiS"^, and hence 

(A.9) P,il-Ar + co^z^)P^<S^P, 

V ' 

all terms positive 



Since [Pi(l - A^)Pi, S^Pi] = and Pi(l - A^)Pi < S^Pi (from (|A9])), we have by Lemma 

Agg 



(A.IO) Pi(l-A.)2Pi<5^Pi 

Since [Pi{-d^, + io^z^)Pi, S^Pi] = and Pi{-d^, + io^z^)Pi < S^Pi (from (|A9|), we have by 
Lemma A.3| 3 ) 

(A.ll) P^{-^l + u'zyP,<S'P^ 

Expanding and "integrating by parts" 

(A. 12) {-dl + u^z^ f = at - 2u^d,z^d, + u^z* +B + B* 

^ V ' 

terms all positive 
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where B =^ —2u'^dzZ. We claim 

(A.13) Pi{B + B*)P,<S'Pi 

Since Wd^PifW^ < \\SPJ\\l2 and uj\\zPJ\\l2 < \\SPJ\\l2, it follows by Cauchy-Schwarz 
that 

,21 T3„/a D f ~D /\l ^ , ,11 CD /l|2 ^ II C2 D ^l|2 



''\Re{d^PJ,zPJ)\<u\\SPJ\\% < \\S'P,jnLi 



which is equivalent to (A.13). By (A. 11), (A. 12), (A.13), we obtain 
(A.14) P^{dt)Pi<S'Pi 
Now, ( |Al0| , dAlil ) imply 

(A.15) Pi(l-A,)'Fi<5'Pi 



Then ( |A.5[ ), ( |A.6[ ), ( |A.7| ), ( |A.8[ ) follow from Sobolev embedding and ( |A.9[ ), ( |A.15[ ). For 
example, to prove (A. 7), we apply 3D Sobolev embedding and ( |A.15 ) to obtain 

l|v,Pi/|U6 < WKPjy. < ws'pju^ < ws'fh.. 

Next we prove 



(A.16) 
(A. 17) 
(A. 18) 
(A.19) 
Recall that 
(A.20) 

We have 
(A.21) 

By Cauchy-Schwarz 
(A.22) 
Also, 
(A.23) 



\\VrPof\\^<^'^'\\Sfh. 

\\Pof\\Le<U:'/'\\Sfh. 

I|V.Po/||l«<u;^/'||5V||lp 

\\Pof\\L^<Oo'/'\\S'fU. 

Pof{x, z) = I f{x, z')K{z') dz' h^{z) = (/(x, ■), h^)^: K{z) 



V,Pof{x,z) = (V,/(x,-),U h^) 

I|v.Po/IIlp < liv./iu. < \\sfh. 



d^Pof{x,z) = {f{x,-),huj) d^K{z) 
and hence by Cauchy-Schwarz, 

(A.24) \\d.P,f\\Ll<uj'"\\f\\Ll 



(A.22) and (A.24) together imply (A.16). By Cauchy-Schwarz, Minkowski, and 2D Sobolev, 

\\Pof\W<W{x,-),K)ALi\\KU 

\\\hui\\Ll\\huj\\Ll 
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Since (1 — A^.) < 5*^, we obtain (A. 17) as a consequence of the previous estimate. Next, we 
prove ( A.18[ ). By (A.21), Cauchy-Schwarz, Minkowski, and 2D Sobolev, 

I|v.Po/IIl«< ll(v./(x,-),/i.)||L|l|/iJU6 

< l|V./|U6^2 a;V6 

< IK-A.P/IUp c.^' 
Since (-A^)^/^ < (1 - AJ^ < ^^^^^^^ 

(A.25) llV.Po/lUe <u;'/15V||l? 



By (A. 23), Cauchy-Schwarz, Minkowski, and 2D Sobolev, 

WdzPofhe < \\{f{x,-),huj)\\L<l\\dzhu,\\Lt 

< ll/llLELi 

< ||(l-A.)i/3/IU? 
Since (1 - A^f/'^ < (1 - A^)^ < 5^ we obtain 

(A.26) \\d.PofU<\\S'f\\^co'/' 



Combining ( |A.25[ ) and ( |A.26[ ), we obtain ( |A.18[ ). Next, we prove ( |A.19[ ). By ( |A.20[ ) and 2D 
Sobolev, 



I^o/IIl- < ||(/(a;,-),^a;)||L-||/i.^||Li 



< ||(1-A.)H7||^, 



Since (1 — Aa;)^+^^ < (1 — A^-)^ < S'^, we obtain (A. 19) as a consequence of the previous 
estimate. 

□ 



Note that combining (|A.5[)-(|A.8[) and ( |A.16[ )-( |A.19[ ) yeilds (|A.l|)-(|A.4[). 

Let 



S=(l-A,. + u;(-l-9,2 + ^2))i/2 

S'>1- Ar 



S'Pi > Pi(l - A^. - tudi + Loz')Pi 



Lemma A. 5. 

(A.27) 
(A.28) 

(A.29) S'Pi > ioPi 

Proof. Directly from the definition of S, we have 

(A.30) Pi(l - A^ - udl + uz^)Pi < ooPi + S^Pi 

^ V ' 

all terms positive 

By spectral considerations 

(A.31) 2cc;Pi < uj{-l -dl + z^)Pi < S'^Pi 
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Combining (A.30) and (A.31) yields (A.28). Also, (A.29) follows from (A.31). Next, we 



establish (A. 27) using (A.28). It is immediate that 
(A.32) > (1 - A J 

On the other hand, since Pq is just projection onto the smooth function e" 



(A.33) 



By (A.28), 



(A.34) 



Pi{-dl)Pi < S^Pi < 



By Lemma |A^3]([T]), ( |X33| ), ( |AlM| ) 
(A.35) 



The claimed inequality (A. 27) follows from (A.32) and (A.35). 
Lemma A. 6. Suppose a : L^(]R^^) — )■ L^(]R^'^) has kernel 



□ 



for some G L^{R^^), and let A, B : L'^{R^'') -> L^{R^''). Then the composition AaB has 
kernel 



{AaB){rk,r'i.) = j {Aip){vk,rN-k){B*'4>){r'f^,rN-k) dvN-k 



It follows that 



TiAaB = {A^,B*ip). 



Let /Cfc denote the class of compact operators on L^(]R^'^), Cl denote the trace class oper- 
ators on L'^{M.^^), and £| denote the Hilbert-Schmidt operators on L'^{M.^^). We have 

^k'^-^k'^ 

For an operator J on L^(]R^^), let \ J\ = {J* J)^^"^ and denote by J(rfc, r'^,) the kernel of J and 
|J|(rfc,r'^) the kernel of |J|, which satisfies |J|(rA;,r'^) > 0. Let 

be the eigenvalues of \ J\ repeated according to multiplicity (the singular values of J). Then 



\J\ 



I op 



op 



\\J\\cl = hnhl = ll^(rfc,r'fc)|U2(r,,r'J = (Tr JV)^/^ 

\\J\\cl = Wf^nhl = lll^l(rfe,rfe)||ii(r,) = Tr |J| 

The topology on ICk coincides with the operator topology, and JCk is a closed subspace of the 
space of bounded operators on L^(]R^'^). 
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Lemma A. 7. Let x be a smooth function on such that xiO — ^ f^^ 1^1 — 1 ^'^^ x(0 — 
for \i\ > 2. Let 

{QMf){v,)= I e'''-^-\[x{M-\^)f\i,)di, 

i=i 

With respect to the spectral decomposition of L^(]R) corresponding to the operator Hj = 
—d1. + let he the orthogonal projection onto the sum of the first M eigenspaces (in 
the Zj variable only). Let 

k 

Rm = W Z^M 

i=i 

(1) Suppose that J is a compact operator. Then Jm == RmQmJQmRm J in the 
operator norm. 

(2) HjJm, JmHj, Ar -Jm and JpiArj are all bounded. 

(3) There exists a countable dense subset {Tj} of the closed unit ball in the space of 
bounded operators on L^{R^'') such that each Tt is compact and in fact for each i 
there exists M (depending on i) such that % = RmQ mTiQ mRm ■ 

Proof. (1) \i Sn ^ S strongly and J G /Cfc, then Sn-J — > SJ in the operator norm and 
JSn — >■ ^5" in the operator norm. (2) is straightforward. For (3), start with a subset {Yn} 
of the closed unit ball in the space of bounded operators on L'^(M?^) such that each is 
compact. Then let {Tj} be an enumeration of the set RmQmYuQmRm where M ranges over 
the dyadic integers. By (1) this collection will still be dense. □ 

Appendix B. Deducing Theorem 11.11 from Theorem 11.21 

The argument presented here which deduces Theorem 1A_ from Theorem [1.2 has been used 
in all the nD to nD work. We refer the readers to them for more details. We first give the 
following proposition. 



1.1 



Let X e 



Proposition B.l. Assume ip^^iO) satisfies (a), (b) and (c) in Theorem 

(M) be a cut-off such that Q ^ x ^ ^, x{s) = ^ for Q ^ s ^ I and x{s) ^ for s ^ 2. 
For K, > 0, we define an approximation ofipj^^lO) by 



This approximation has the following properties: 
(i) V'Artj(O) verifies the energy condition 



(u) 



sup 

N,uj 



L2 
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(Hi) For small enough k > 0, "^tvi^IO) ^■^ asymptotically factorized as well 



lim Tr 

7V,tj— s>oo 



7^^^^(0,a;i,2;i;xi,2;i) - (f)Q{xi)(f)Q{x[)h{zi)h{z[) 



0, 



where 7^^^^^ (0) is the marginal density associated with "ipN^iO), and 0o ihe same as in 
assumption (h) in Theorem \l.l\ 



Proof. Proposition B.l follows the same proof as [26, Proposition 9.1] if one replaces by 
[Hjsi^u] — Ncu) and Hn by 

^ 1 



J=2 



l<i<j<N 



Via (i) and (iii) of Proposition 1.2, ipj^^^iO) verifies the hypothesis of Theorem 



enough k > 0. Therefore, for (t) , the marginal density associated with e^^^^^^^ipN ui^) 
Theorem 1.2 gives the convergence 



1.2 



□ 



for small 



(B.i: 



lim Tr 

AfjOj— s>oo 



7^,^^(^,Xfc,Zfc;x'fc,z'fc) -Y[(j){t,Xj)(j){t,x'j)hi{zj)hi{zj] 



for all small enough k, > 0, all k ^ 1, and all t G M. 

For 7^1^ (t) in Theorem 1.1, we notice that, VJ^''^ G JCk, Vt G M, we have 



I + 11. 



Convergence (B.l) then takes care of II. To handle I , part (ii) of Proposition 1.2 yields 



L2 



L2 



which implies 



Tr/'^) (7!J,lW-7^2^W)|^C7||j«||^^,|. 



Since k > is arbitrary, we deduce that 



lim 

N,uj—>oo 



0. 



i.e. as trace class operators 

7^,1 it) ^ 10 it) ® h) (0 (t) ® /iil^'^ weak*. 
Then again, the Griimm's convergence theorem upgrades the above weak* convergence to 



strong. Thence, we have concluded Theorem 1.1 via Theorem 1.2 and Proposition B.l 
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